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Abstract. 

We develop the reconstruction program for the number of modified gravities: scalar- 
tensor theory, f{R), F{G) and string-inspired, scalar-Gauss-Bonnet gravity. The 
known (classical) universe expansion history is used for the explicit and successful 
reconstruction of some versions (of special form or with specific potentials) from all 
above modified gravities. It is demonstrated that cosmological sequence of matter 
dominance, decceleration-acceleration transition and acceleration era may always 
emerge as cosmological solutions of such theory. Moreover, the late-time dark energy 
FRW universe may have the approximate or exact ACDM form consistent with three 
years WMAP data. The principal possibility to extend this reconstruction scheme to 
include the radiation dominated era and inflation is briefly mentioned. Finally, it is 
indicated how even modified gravity which does not describe the matter-dominated 
epoch may have such a solution before acceleration era at the price of the introduction 
of compensating dark energy. 
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1. Introduction 

The explanation of tlie current universe speed-up which is due to mysterious dark 
energy remains to be one of the most fundamental problems of modern cosmology 
and theoretical physics. The number of theoretical models have been developed 
aiming to describe the dark energy universe (for recent review, see [H [2]). The 
obtained observational data indicate that current universe has the effective equation 
of state parameter w being very close to —1 (for the review of observational data from 
the theoretical point of view, see and for description of observable cosmological 
parameters, see ^). Hence, it is still not quite clear if the universe lives in the phantom 
era {w less than —1), in the (most probably) cosmological constant epoch {w = — 1) 
or in the quintessence phase {w more than —1). In such a situation, even phantom 
cosmological models (see [S, El [7] and refs. therein) which show quite strange properties 
are not ruled out because even being not phantomic one, the universe may currently 
enter to phantom dominated epoch. 

The dark energy problem may have the gravitational origin, calling for possible 
late-time modification of General Relativity. The number of modified gravity models 
(for a review, see [8j) starting from the simplest 1/R theory [9l[l0] have been proposed 
as gravitational alternative for dark energy. Such models may be inspired by string/M- 
theory considerations [HI |12]. They may lead to (phantom, cosmological constant 
or quintessence) effective cosmology at late times without the need to introduce the 
scalar fields with strange properties (like negative kinetic energy) or with complicated 
potentials. 

Generally speaking, modified gravity looks very attractive as it gives the qualitative 
answers to the number of fundamental questions about dark energy. Indeed, the origin 
of dark energy may be explained by some sub-leading gravitational terms which become 
relevant with the decrease of the curvature (at late times). Moreover, there are many 
proposals to consider the gravitational terms relevant at high curvature (perhaps, due to 
quantum gravity effects) as the source of the early-time inflation. Hence, there appears 
the possibility to unify and to explain both: the inflation and late-time acceleration 
as the modified gravity effects. Meanwhile, at the intermediate epoch the gravity may 
be approximated by General Relativity. Similarly, the coincidence problem may be 
explained. It is remarkable that modified gravity may pretend to describe also dark 
matter. 

There are various ways to describe the modified gravity theory. Not only f{R) 
but also modified Gauss-Bonnet theory [13] may be presented in the form of scalar- 
tensor gravity (perhaps, with higher derivative terms). Note that scalar-tensor theory 
itself may be understood as some kind of modified gravity. Moreover, one can present 
modified gravity as General Relativity with quite complicated ideal fluid having 
non-standard equation of state. As it has been explained in ref. [15], this mathematical 
equivalence of different representations does not mean the physical equivalence (due to 
transformed form of the physical metric tensor). Nevertheless, it is remarkable that 
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modified f{R) gravity has some close analogy with non-trivial f{p) equation of state 
ideal fluid [HI [7]. In general, modified gravity may be considered as General Relativity 
(GR) with inhomogeneous equation of state ideal fluid [17]. (The particular example of 
this sort is given by ideal fluid with time- dependent bulk viscosity jT8]). 

Of course, modified gravity should be fitted against the observational data (as is 
done in refs. |19l [20] , mainly, for f{R) theory). Similarly, it should pass the Solar System 
tests (for f{R) models this is widely discussed in refs. [211 1221 |23l |20] while SdS metrics 
are studied in [21]). In principle, it may pass the requested tests as some cancellation 
of the contributions between different terms may occur as in the model of ref . |23j . 
Moreover, the corresponding initial (boundary) conditions of the form of the action at 
some (current) time of the universe expansion history may be used to bring the action 
to the approximate GR form (with possible earlier and (or) future deviations from it). 

In such a situation, any realistic classical gravitation should describe the known 
sequence of the cosmological epochs (perhaps, without the inflationary era where 
quantum gravity effects may be quite strong). The purpose of this work is to 
review several versions of modified gravity and to show how it may be reconstructed 
from the known universe expansion history (for recent discussion of reconstruction 
program for usual gravity with matter, see [25]). In the next section this problem 
is solved for two scalar-tensor theory. The reconstruction method is developed. 
Several examples of realistic cosmology are considered: oscillating universe, ACDM 
cosmology and asymptotically de Sitter space. The corresponding scalar potentials 
admitting the sequence of matter dominated phase, decceleration-acceleration transition 
and acceleration are constructed in all above cases. Section three is devoted to the 
reconstruction of f{R) gravity following ref. [15l |26]. The modified gravity admitting 
the unification of matter dominated and acceleration era is formulated. Its asymptotic 
behaviour at early and late times is found. It is shown that it does not conflict with 
three years WMAP data. 

In section four the alternative scenario for the model of ref. [23] is developed. Adding 
the compensating dark energy which role is negligible at the acceleration epoch it is 
shown that modified gravity [23] may be viable (the matter dominated phase occurs 
before the cosmic acceleration). In section five the reconstruction of string- inspired 
scalar-Gauss-Bonnet gravity and of modified Gauss-Bonnet gravity is made. Such 
theories have been suggested as dark energy models in refs. [T2| [13], respectively. Its 
versions where dark energy universe emerges after matter dominance are given. Some 
discussion and outlook are presented in the Discussion section. 
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2. The reconstruction of scalar-tensor theory from the expansion history of 
the universe 

2.1. General formulation 

In this subsection, we show how any cosmology can be reproduced in the scalar-tensor 
theory (or scalar-tensor theory may be reconstructed). Our consideration is based on 
the method developed in ref. [27] (for recent study of scalar-tensor cosmology applied 
for dark energy description see [281 E] and for earlier attempts, see |29]). Note that 
studying the transition from non-phantom phase to phantom phase, the instability of 
the one scalar model becomes infinite at the transition point [27j. That is why we only 
consider two scalar-tensor theory, whose instability could be always finite. Moreover, 
with two scalars it is easier to fit with universe history expansion. 
We now consider two scalar-tensor gravity 

S = ld''xy^!^^R-^u;i<P)d,<pd'^<P-^vix)d,xd''x-V{<P,x)} • (1) 

Here uj{(f)) is a function of the scalar field (p and r]{x) is a function of the another scalar 
field X- It is assumed the spatially-fiat FRW metric ds"^ = —dt^ + a{tY Yh=i {dx'') and 
that and x only depend on the time coordinate t. The FRW equations lead to 

a;(0)02 + rj{x)x' = -^H , \/(0, x) = ^ i^H^ + h) . (2) 

Then if 

^it) + nit) = -^J'{t) , \/(t, t) = 1 (3/(t)2 + fit)) , (3) 

the explicit solution follows 

= X = t, H = f{t). (4) 

One may choose that u should be always positive and rj be always negative, for example 

2 



Vix) = - ^V'a2 + f(x)2<0. (5) 
Here o; is a constant. Define a new function /(^, x) by 

m X) = -y (/ #^(0) + / dxvix)) ■ (6) 
The constant of the integration could be fixed to require 

f{t,t) = f{t). (7) 
If Vi<P, X) is given by using /(0, x) as 

V^(0, X) = -, (3/(0, X) + + 1 , (8) 
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the FRW and the scalar field equations are also satisfied: 







+ 3Huj 



+ 



= V{X)X + lv'{x)x' + ^Hr^{x)x + ^^^^ 
2 dx 



(9) 



We now investigate the (in)stability of the model. By introducing the new 
quantities, X^, X^, and Y as 

mx) 



X 



Y 



the FRW equations and the scalar field equations 
dX^ u'{(f)) ( 2 



(10) 



are: 



dN 
dX^ 



2Huj{(f)) 

v'ix) 

1 



dN 
dY 

dN ~ 2^FlP 
Here d/dN = H^^d/dt. In the solution 



(X|-1)-3(X^-F) 



{X4x<^r-i) + x^(x^r-i)} 

), Xd, 



perturbation may be considered 
X^ = 1 + 6Xj) , 



X^ = 1 + 5X^ , 



Y 



X, = y 



6Y . 



(11) 

1. The following 
(12) 



Hence 



d 

dN 



M 



( sx^ 


] 




( sx^ 








= M 






\ 6Y 


) 




6Y 


) 



Hijj{(t>) 


1 



-3 



V(x) 

Hri(x) 
1 



-3 



3 
3 



(13) 



The eigenvalues of the matrix M are given by solving the following eigenvalue equation 







Hvix) 



X 



1 



2K2/f2 



A + 



v'ix) 

Hvix) 



(14) 



2k'^H^ \' ' Huji(j)) 

The eigenvalues (fH|) for the two scalar model are clearly finite. Hence, the instability 
could be finite. In fact, right on the transition point where H = f'{t) = and therefore 
/'(0) = fix) = 0, for the choice in ([5]), we find 



2a 

^(0) = -Vix) = — , ^'(0) 



2H 



v'ix) = . 



(15) 



Then the eigenvalue equation flMj) reduces to 



= + i-A -B + 6)\^ + iAB - 3A - 35 + 9) A - -AB + 95 



A = 



H 



B = 



a 



(16) 
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Here we have chosen a > 0. Then the eigenvalues are surely finite, which shows that even 
if the solution (jl]) is not stable, the solution has non-vanishing measure and therefore 
the transition from non-phantom phase to phantom one can occur. We should also note 
that the solution (jl]) can be in fact stable. For example, we consider the case A,B ^ 0. 
Eq. (|T6l) further reduces to 

= A(A + 3)^ . (17) 

Then the eigenvalues are given by and —3. Since there is no positive eigenvalue, the 
solution (jll) is stable in the case. 

It is not difficult to extend the above formulation to the multi-scalars model, whose 
action is given by 



:i8) 



Here uJi{(f)i) is a function of the scalar field (pi. We now choose Ui to satisfy 

= -A/'W (19) 

I 

by a proper function, and we also choose the potential Vi^cpi) by 
Here 

2 

/(0,)^-^5:|rf0,cu,(0,). (21) 
The constant of the integration in fl2Tl) is determined to satisfy 

Then a solution of the FRW equations and the scalar field equations are given by 

0, = t , H = fit) . (23) 
The rest consideration coincides with the one given for two-scalars model. 

2.2. The examples of scalar-tensor theory reconstruction 

We now consider some explicit examples in order to demonstrate how the scalar-tensor 
theory can be reconstructed from known expansion history of the universe using the 
above developed formulation. 

As a first example, we consider the model given by 

/(0) = h^ + hi sin(z/0) . (24) 
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Here it is assumed ho, hi, and u are positive. By choosing a = hiv in one finds 



uj{(t)) = |cos(z/0) - y^l + cos2(z/(/ 

^(X) = -^K^\/lTcOs2(^ , 

X) = /io + h sin(z/0) - ^ {i^ (l/v^, ^^0) - E {1/V2, ux) } . (25) 
Here E{k,x) is the second kind eUiptic integral defined by 



E{k,x)=l da;Vl - Psin^x . (26) 
Jo 

Note that similar reconstruction in case of the single scalar-tensor theory was presented 
in [27]. Eq. flMl) shows that the Hubble rate H is given by 

H = ho + hi sin(z/t) , (27) 

which is oscillating. When h^ > hi > 0, H is always positive and the universe is 
expanding. Since 

H = hiucos{ut) , (28) 
when hiu > 0, Wes, which is defined by 

= r - 1^ • '''' 

is greater than —1 (non-phantom phase) when 

2n - -) TT < z/t < f 2r2 + vr , (30) 



2/ V 2 

and less than —1 (phantom phase) when 



2n + -]'K<ut< (2n + -]7i. (31) 



2/ V 2. 

In f l30p and (1311) . n is an integer. Hence, in the model (1241) . there occur multiply 
oscillations between phantom and non-phantom phases. It could be that our universe 
currently corresponds to late-time acceleration phase in such oscillatory regime (for 
recent reconstruction examples for oscillatory universe, see [30] and references therein). 

The present universe is expanding with acceleration. On the other hand, there 
occured the earlier matter-dominated period, where the scale factor a behaves as 
a ~ t^/^. Such behavior could be generated by dust in the Einstein gravity. The baryons 
are dust and (cold) dark matter could be a dust. The ratio of the baryons and the dark 
matter in the present universe could be 1 : 5 or 1 : 6, which should not be changed even 
in matter dominant era. It is not clear what the dark matter is. For instance, the dark 
matter might not be the real matter but some (effective) artifact which appears in the 
modified/scalar-tensor gravity. In the present section, it is assumed that not only dark 
energy but also the dark matter may originate from the scalar field 0. 

We now investigate that the transition from the matter dominant period to the 
acceleration period could be realized in the present formulation (for single scalar-tensor 



Modified gravity and its reconstruction from the universe expansion history 



8 



theory, see [21] )• In the following, the contribution from matter could be neglected since 
the ratio of the matter with the (effective) dark matter could be small. 
First example is 

H = fit) = g\t) =90 + J. (32) 

When t is large, the first term in fl5^ dominates and the Hubble rate H becomes 
a constant. Therefore, the universe is asymptotically deSitter space, which is an 
accelerating universe (for recent examples of late-time accelerating cosmology in scalar- 
tensor theory, see E]). On the other hand, when t is small, the second term in (15^ 
dominates and the scale factor behaves as a ~ t^^ . Therefore if gi = 2/3, the matter- 
dominated period could be realized. Here, one of scalars may be considered as usual 
matter. 



Since 



fit) = -y < , (33) 



instead of ([5]), by using a positive constant a, we may choose 



u;{<j>) = -^ii±^/'(0) > , vix) = ^/'(0) < , (34) 



that is 



^(<^) = — 7^72' ^W = 7^7I- (35) 



2( 1 + a) gi ^^^^ ^ 2a 

One should note that in the limit a — 0, the scalar field x decouples and single scalar 
model follows. Then one obtains 

'l + a)gi agi 



/(0,X) = ^0 + 



X 



Hence, for scalar-tensor theory with such potentials the matter dominated era occurs 
before the acceleration epoch. 

Before going to the second example, we consider the Einstein gravity with 
cosmological constant and with matter characterized by the EOS parameter w. FRW 
equation has the following form: 

Here P is the inverse cosmological constant. The solution of fl37j) is given by 

a = aoe^(*) , g{t) = ^^^^ In (^a sinh (^^^^ (t - to)) ) • (38) 
Here to is a constant of the integration and 

2 2/2 -3(l+w) /onN 

a =-K I Potto . (39) 
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As the second example, we consider metric 



m = g'it) - 
fit) = g"{t) 



1 



coth 



'3(1 +«;) 



21 



3(l + w) 



2/2 

it is convenient to use (!34l) instead of (l5| 
.(^) ^3(1 + + a) 



sinh 



for scalar-tensor theory. Since 
it-t^)^ , 



21 



it 



< 



(40) 



Then one arrives at 

'm + w) 




(41) 



Thus, in both examples, (!32|) and (!38|) . there occurs the matter dominated stage, the 
transition from matter dominated phase to acceleration phase and acceleration epoch. 
In the acceleration phase, in the above examples, the universe asymptotically approaches 
to deSitter space. This does not conflict with WMAP data. Indeed, three years WMAP 
data have been analyzed in Ref. [32]. The combined analysis of WMAP with supernova 
Legacy survey (SNLS) constrains the dark energy equation of state wde pushing it 
towards the cosmological constant. The marginalized best flt values of the equation of 
state parameter at 68% confldance level are given by —1.14 < wde < —0.93. In case 
of a prior that universe is flat, the combined data gives —1.06 < wde ^ —0.90. In 
the examples ( l32i) and ( 138|1 . the universe goes to asymptotically deSitter space, which 
gives wde —I, which does not, of course, conflict with the above constraints. Note, 
however, one needs to flne-tune go in (l32l) and 1// in ( l38l) to be (70 ~ 1/^ ~ 10~^^ eV, in 
order to reproduce the observed Hubble rate Hq ~ 70 kms^^Mpc^^ ~ 10^^'^ eV. 

The flnal remark is in order. There is no problem to include into the above 
reconstruction scenario usual matter (say, ideal fluid, dust, etc) and (or) to flnd the 
scalar potentials corresponding to realistic cosmology in multi-scalar case. Moreover, in 
the same way one may include the radiation dominated epoch where quantum effects 
may be still neglected in the above scenario. In the next section it is shown how similar 
reconstruction scheme may be developed for modifled gravity. 
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3. Reconstruction of modified gravity with the unification of matter 
dominated and accelerated phases 

3.1. General formulation 

In the present section we develop the general reconstruction scheme for modified gravity 
with f{R) action. It is shown how any cosmology may define the implicit form of the 
function /. The starting action of modified gravity is: 

S = J d''x{f{R)+C^.,tter} . (42) 

First we consider the proper Hubble rate H, which describes the evolution of the 
universe with matter dominance era and accelerating expansion (for discussion of various 
accelerating cosmologies from above action, see [231 [331 [9l [10]). It turns out that one can 
find /(-R)-theory realizing such a cosmology (with or without matter). The construction 
is not explicit and it is necessary to solve the second order differential equation and 
algebraic equation. It shows, however, that, at least, in principle, we could obtain any 
cosmology by properly reconstructing a function f{R) on theoretical level. 
The equivalent form of above action is 

S = J d'xy^ {P{<P)R + Q{<f)) + . (43) 

Here P and Q are proper functions of the scalar field and ^matter is the matter 
Lagrangian density. Since the scalar field does not have a kinetic term, it may be 
regarded as an auxiliary field. In fact, by the variation over 0, it follows 

O = P'(0)i? + g'(0) , (44) 
which may be solved with respect to 0: 

= 0(i?) . (45) 
By substituting (H5l) into (H3ll . one obtains /(-R)-gravity: 

S = J d\^ {f{R) + C^^tter} , f{R) = P{<I>{R))R + Q{<P{R)) . (46) 
By the variation of the action fj43|) with respect to the metric (7^1,, we obtain 

= -^g,, {P{<P)R + Q(0)}-i?^,P(0)+V^V,P(0)-^7,.V2p(0) + ^T^, .(47) 
The equations corresponding to standard spatially-fiat FRW universe are 

= - QH'P{<P) - Q(0) - QH^^^^T^ + p , (48) 

dt 

= (4H + 6H^) + Q(^) + + W^-^MA + p . (49) 

By combining (HTj) and (148|) and deleting Q{(f)), we find the following equation 

= 2^^-2h'-^^4HPW^,^,. (50) 

dt'^ do 
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As one can redefine the scalar field properly, we may choose 

= t. (51) 

It is assumed that p and p are the sum from the contribution of the matters with a 
constant equation of state parameters Wi. Especially, when it is assumed a combination 
of the radiation and dust, one gets the standard expression 

p = proa + pdoa , p = —a , (52) 

with constants pro and pdo- If the scale factor a is given by a proper function g{t) as 

a = aoe^^*) , (53) 
with a constant oq, Eg. (1491) reduces to the second rank differential equation: 



+ 5: (1 + w,) p,oao '^'"'"'^e-3(^+'"»)^(^) . (54) 

i 

In principle, by solving fl54p we find the form of P{(f)). Using (HSj) (or equivalently dH])), 
we also find the form of Q{4>) as 

g(0) = -6(^7'(0))^P(0)-6^'(0)^ + 5:Aoao'^'-""^^e-=^(^+-)^W .(55) 

Hence, in principle, any cosmology expressed as flS5]) can be realized by some specific 
/(i?)-gravity. 

3.2. f{R) gravity: the transition of matter dominated phase to the acceleration phase 

Let us consider realistic example where the total action contains also usual matter. The 
starting form of g{(j)) is 

^7(0) = M0) In (^|-^ , (56) 

with a constant 0o- It is assumed that /i(0) is a slowly changing function of 0. We use 
adiabatic approximation and neglect the derivatives of h{4>) {h\(j)) ~ /i"(0) ~ 0). Then 
the solution of Eq. (15^ has the following form [151 EE] : 

P(0) = + W + ^ pi(0)0-3(i+-«)M0)+2 _ (57) 

i 

Here p± are arbitrary constants and 



/i(0)-l± J/l(0)2 + 6/1(0) + 1 

n±{.(p) = '-^ , 

pm - -{il + w)p,oao'''^-^UT^'"'''''} 

X {6(1 + w){4: + 3w)h{(l)f - 2 (13 + 9w) h{<p) + 4}"^ . (58) 
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Especially for the radiation and dust, one has 

~ 3a^ (4O/z(0)2 - 32/i(0) + 4) ' 

^'^'^^ = " al (24/i(0)2-26/i(0)+4) ' ^^^^ 
We also find the form of Q{4>) as 

Q(0) = - 6/i(0K m) + ^+(0)) 0"^^^^"' 
- 6h{(f))p_ {h{(f)) + n_(0)) 
+ ^ {-6/i(0) (-(2 + 3w)/i(0) + 2) 

+P.oao '^'+"Vf • (60) 

Eq.dSnD tells that H ~ h{t)/t and i? ~ 6 i-h{t) + 2h{ty) /t^. Let assume lim<^_,o ^^(0) = 
hi and lim^^oo ^(0) = hj. Then if < /ij < 1, the early universe is in decceleration 
phase and if hf > 1, the late universe is in acceleration phase. We may consider the 
case ~ hm is almost constant when (p ^ tm (0 <C ■'^ +oo). If hi, hj > 1 and 
< hm < the early universe is also accelerating, which could correspond to the 
inflation. After that the universe deccelerates, which corresponds to matter-dominated 
phase with h{(j)) ~ 2/3 there. Furthermore, after that, the universe could be in the 
acceleration phase. 

The simplest example is 

hi + hfq(f 

M0) = , (61) 

with constants hi, hf, and q. When — 0, /i(0) hi and when — » oo, /i(0) — >■ hf. If 
q is small enough, /i(0) can be a slowly varying function of 0. Then we find [26] 

02 = $o(i?) , ^o^aV' + ai^' , 
a± = 



/9o = 



-/^o±\//^g-y 
2 

2 (2R + 6hfq - 12hjq 



27g3/?3 

(2/? + 6/i/g - 12h}q) {R + G/z^g + 6/i/g - 4/ii/i/g) 

3g^ 

+ 6hi - 12h^i , 

{2R + 6hfq-12h}qy R + 6h,q + 6hfq - Ahihfq , , 

" ^ • ^''^ 

There are two branches besides $o but the asymptotic behaviour of R indicates that 
we should choose $o in ( |62l) . Then explicit form of f{R) could be given by using the 
expressions of -P(0) (!57|) and (5(0) (!60l) as 

/(i?) = P (V<J^o(i?)) + g (\/'^'o(i?)) . (63) 
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One may check the asymptotic behavior of f{R) in For simphcity, it is 

considered the case that the matter is only dust {w = 0) and that p_ = 0. Then 

P(0) = + . (64) 

One may always get — {—3h + 2) > in Here n+ is defined in (1581) . Then when 
is large, the first term in flM|) dominates and when is small, the last term dominates. 
When (f) is large, curvature is small and 0^ ~ 6 {—hf + 2hf) / R and h{(f)) h{oo) = hf. 
Hence, Eq. fl6^ shows that 

P(0)~p+( ^ '-j , (65) 

and therefore f{R) ~ ^-(M</')-5+v/^J+6X7TT)/4_ Especially when /i > 1, we find 
/(i?) ~ i?-'^//^^ Therefore there appears the negative power of R (for review of such 
theories, see [8]). As H hf /t, if hf > 1, the universe is in acceleration phase. 

On the other hand, when curvature is large, we find 0^ ~ 6 {—hi + 2hi) / R and 
h{(j)) h{0) = hi. Then (!57|) shows P{(f)) ~ pd(O)0"^'^*+^. If the universe era 
corresponds to matter dominated phase {hi = 2/3), P{4>) becomes a constant and 
therefore f{R) ~ R, which reproduces the Einstein gravity. 

Thus, in the above model, matter dominated phase evolves into acceleration phase 
and f{R) behaves as f{R) ~ R initially while f{R) ~ ji-{hW-5+^h-j+6hf+i)/A 
late time. Moreover, for some parameter choice the asymptotic limit of above theory 
reproduces the model of ref . [23] . 

In our model, we can identify 

2,2 

wde = -1 + ^, or hf = —— -, (66) 

3hf 3(1+m;de) 

which tells hf should be large if hf is positive. For example, if wbe = —0.93, 

hf ~ 9.51 ■ ■ ■ and if wde = —0.90, /i/ ~ 6.67 ■ ■ ■. Thus, we presented the example of f{R) 

gravity which describes the matter dominated stage, transition from decceleration to 

acceleration and acceleration epoch which is consistent with three years WMAP. Adding 

the radiation permits to realize the radiation dominated era before above cosmological 

sequence as is shown in [26]. 



4. Modified gravity and compensating dark energy 

In this section we present another approach [15] to realistic cosmology in modified 
gravity. Specifically, we discuss the modified gravity which successfully describes the 
acceleration epoch but may be not viable in matter dominated stage. In this case, 
it is demonstrated that one can introduce the compensating dark energy (some ideal 
fiuid) which helps to realize matter dominated and decceleration-acceleration transition 
phases. The role of such compensating dark energy is negligible in the acceleration 
epoch. 
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We now start with general /(i?)-gravity action fl42p . In the FRW metric with flat 
spatial dimensions one gets 

p = f(R)-6{H + H'-H^\ f\R) , 




p = - f{R) -2[-H-3H' + ^^+ 2Hj^ f{R) . (67) 

Here R = 6H + 12H^. If the Hubble rate is given (say, by observational data) as 
a function of t: H = H{t), by substituting such expression into fl^7|) . we find the t- 
dependence of p and p as p = pit) and p = pit). If one can solve the first equation with 
respect to t as t = t(p), by substituting it into the second equation, an equation of state 
(EOS) follows p = p (tip)). Of course, p and p could be a sum with the contribution of 
several kinds of fluids with simple EOS. 

We now concentrate on the case that f{R) is given by [23] 

/(^) = -|r + ^ + ^^'- (68) 

Furthermore, we write H{t) as H{t) = h{t)/t and assume h{t) is slowly varying function 
of t and neglect the derivatives of h{t) with respect to t. 

First we consider the case that the last term in ( l68l) dominates f{R) ~ which 
may correspond to the early (inflationary) epoch of the universe. It is not difficult to 
find 

3Q/3 i-l + 2h{t)) hjt)^ 36(3{-l + 2h{t))h{3h{t) + l) 

If h goes to infinity, which corresponds to deSitter universe, we find p ^ p ^ h^ although 
R ^ h^. Therefore p, p <^ jSR'^ and contribution form the matter could be neglected. 
Then the infiation could be generated only by the contribution from the higher curvature 
term. 

Second, we consider the case that the second term in (168!) dominates /(-R) ~ R/2k^, 
which may correspond to the matter dominated epoch after the infiation. In this case 
p and p behave as 

I2h{t) + 6/i(t)2 Ahit) - 6h{t)^ 

P ^2 ' P ~ 272 • ^'^> 

In the matter dominated epoch, we expect h ~ 2/3 (^a~ti^. Hence, one gets 
p ~ 32/3fi;^t^, p ~ 0. Therefore in the matter sector, dust with = (p = 0) should 
dominate, as usually expected. 

Finally we consider the case that the first term in (l68l) dominates f{R) ~ —a/R"", 
which might describe the acceleration of the present universe. The behavior of p and p 
is given by 

p ~ a {6 (n + 1) (2n + 1) h{t) + 6 (n - 2) h{tf] 
X {-6/i(t) + 12/i(t)2}"""^t2n ^ 

p ~ a [-An [n + 1) {2n + 1) - 2 (Sn^ + 5n + 3) h{t) - 6 (n - 2) /i(t)^} 
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X {-6h{t) + 12/i(t)2}^"~^t2n _ ^7^^) 

Thus, the effective EOS parameter wi = p/p is given by 

wi ~ {-4n {n + 1) {2n + 1) - 2 (Sn^ + 5n + s) h{t) 

-6 {n - 2) h{ty} {6 (n + 1) {2n + 1) + 6 {n - 2) h{tfy^ . (72) 

In order that the acceleration of the universe could occur, we find h > 1. Let us now 
assume that h{t) hf when t —>■ oo. Then one obtains wi Wf. Here 

Wf = [-An {n + 1) {2n + 1) - 2 (Sn^ + + s) hf 

-6 (n - 2) /ij} {6 (n + 1) (2n + 1) /ij + 6 (n - 2) /ij}"^ , (73) 

and H{t) — > hf/t. Since the matter energy density p^^ with the EoS parameter Wf 
behaves as 

p^^ cx a-3(^+"'^) oc exp (^-3(1 + wj) J dtH{t)^ , (74) 

the energy density is 

p^f oc _ (75) 

Comparing (1751) with ( 17Ti) . we find 2n = —3(1 + Wf)hf, which can be confirmed directly 
from (1731). 

From the above consideration, we find p and p contain mainly contributions from 
dust with w = 0, Pdif), pit) = and "dark energy" with w = wi in (!72|) . pi{t), piit). In 
the expressions of p{t) and p{t) in (|67l) . there might be a remaining part: 

PR{t) = pit) - pdit) - plit) , pnit) = Pit) - piit) , (76) 

which may help the transition from the matter dominated epoch to the acceleration 
epoch. By deleting t in the expression of (l76ll . we obtain the EOS for the remaining 
part: 

PR = PnipR) , (77) 
which may be called the compensating dark energy. More concretely, one may have 

Prf-^e , (78) 

and according to (|7T]) . 

~ a {6 (n + 1) (2n + 1) /ij + 6 (n - 2) h'j} 

X {-6hf + 12/.;}—^ t2n^-3(i+».)/; ^^g^ 

In ( 179|l . we choose ti to be large enough. When t ~ to, pit) ~ pd and when t ^ oo, 
pit) ~ Pi. Thus, pr only dominates after t = ti but it becomes smaller in late times. 
Hence, the role of pr (which perhaps may be identified partially with dark matter) is 
only to connect the matter dominated epoch to the acceleration epoch. The number of 
modified gravities can be made cosmologically viable via such scenario. 
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5. Reconstruction of scalar-Gauss-Bonnet theory and F(G')-gravity 



5.1. Scalar-Gauss-Bonnet gravity 



S = d^x^/^ 



(80) 



In this section we show how string-inspired, scalar- Gauss- Bonnet gravity may be 
reconstructed for any requested cosmology. The starting action is 

- R 1 
.2^ ~ 2* 

Here G is the Gauss-Bonnet invariant and the scalar field is canonical in ( IHOi) . Note 
that scalar may be considered as matter component. We assume the FRW universe and 
the scalar field only depending on t. The FRW equations look like p^: 

Q = -^H' + + + 2AH''^^^^ , (81) 
2 dt 

= 1 (2i/ + 3tf^) + - V(4,) - Sff^^i^ 

''^,,HH'Mm.,,„.'Mm, (82) 

dt dt ^ ' 

and scalar field equation: 

O = + 3/f0 + V'(0)+eU0)G . (83) 
Now G = 24 (hH'^ + H^). Combining ([HI]) and ([H2]), one gets 



dt^ dt dt 

2H + ^^-8al(^^»V (84) 



dt \ a dt j 
Eq. fl84l) can be solved with respect to ^i(0(t)) as 

1 /S. 



W{t) ^ f ^){^2H{ti) + ktif) . (85) 
Combining (IHTj) and (185|) . the scalar potential ^^(^(t)) is: 

V{m) = -,H{tf - h{tf - 3a{t)H{t)W{t) . (86) 

We now identify t with /(0) and with g'it) where / and g are some functions. Let us 
consider the model where V^(0) and ^i(0) may be expressed in terms of two functions / 
and g as 



^JM = - / dcPi^^^^ ^f/(0i) 

8 J a'(f(6.)? ' 



9'{f{<Pi)r 

m ^ f #i/'(0i)e-^(^(^^)) i^' {f{<P,)) + . (87) 
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By choosing V{(j)) and ^i(0) as (157|) . we can easily find the following solution for Eqs. (15Ti) 
and ([82]): 

(j) = r\t) (t = /(</>)), a = aoe^^'^ iH = g'{t)) . (88) 

One can straightforwardly check the solution ( l88l) satisfies the field equation ( l83i) . 
An interesting cosmological example is 

^ (l)''e5o* ^ ^ = j-i(i) = 0„inl . (89) 
\to/ ^0 
Here to, go, gi, and 0o are constants. We further choose 



(90) 



Then [/ is a constant: U = Uq and 



,.„).^r"°*f*o+?ir7^)e-. (91) 



tiJ Vto 



Eq.dlO]) shows 



H = go+^j. (92) 

Hence,, when t is small, the second term in (p2!) dominates and the scale factor behaves 
as a ~ Therefore if (^i = 2/3, the matter-dominated period could be realized. On 
the other hand, when t is large, the first term in fl92p dominates and the Hubble rate H 
becomes a constant. Therefore, the universe is asymptotically deSitter space, which is 
an accelerating universe. As in our model the universe goes to asymptotically deSitter 
space, we find wde —1- Therefore our model can easily accommodate three years 
WMAP values of wde- For example, if go ^ 40, one has wde = —0.98. 

In the limit Uq — ^ 0, the Gauss-Bonnet term in fl80l) vanishes and the action 
reduces into that of the usual scalar tensor theory with potential 

2 

V"" to J kHq 

which reproduces the result ^Ij. 

Let us reconstruct the scalar- Gauss-Bonnet gravity from FRW cosmology (138|) . If 
a function g{t) is given by ( l38ll and /(</>) is given by 



y^^) = -2 [9^ + T'-'^n - 4^e-^^/^o , (93) 



/(0) =to- ^ Intanh , (94) 

U{(p) flHTj) becomes a constant again, U = U{(f)). Then V{(p) and ,^(0) are found to be 
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3(1 + w) . 2 /'^V3(l + w) 



2Pk 



3f/n 



cosh 



Ki/3(1 + w 



1 /kJ3(1 + u; 

X < sinh — 

a I 2 




6(0) 



° / dcpi cosh ^ ( - 

l + w)-' \ 



-2/(3(l+t«)) 



kJ3(1 + w) 




1 /m/3(1 
- — sinh , 
a \ 2 



-2/(3(l+«>))-l 



(95) 



s2/(3(l+t«)) 



which again reproduces the result of ref.|3T] in the hmit of Uq 0. 

Eg. (1381) shows that when t ~ tQ, the scale factor behaves as a ~ (t — tpj 
Therefore if w = 0, the matter-dominated period could be realized. On the other hand, 
when t ^ oo, a behaves as a ~ e*/', which tells the universe goes to asymptotically 
deSitter with w^e ^ — 1. Hence, it could be consistent with WMAP data. The FRW 
cosmology of these and another versions of scalar- Gauss-Bonnet gravity may be studied 
in the same way as it was done in refs. |34[ [35] . The investigation of cosmological 
perturbations can be also done [36] in the above model admitting the cosmological 
sequence of matter dominance, decceleration-acceleration transition and acceleration. 



5.2. F{G) -gravity 

The formulation of the previous section can be extended to so-called F(G)-gravity [13] . 
whose action is given by 

R 



S 



(96) 



The above action could be rewritten by introducing the auxilliary scalar field as 

- R 
2^ 



r(0)-ei(0)G 



(97) 
(98) 



S = J d xy— ^ 
By the variation over (p, one obtains 

which could be solved with respect to as 

= 0(G) . (99) 

By substituting the expression (199!) into the action (l97j) . we obtain the action of F{G)- 
gravity with 

F{G) = -V (0(G)) + 6 (0(G)) G . (100) 
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Note that the action (P7|) could be obtained also by dropping the kinetic term of in 
the action (|80D . 

Assuming the spatially-flat FRW universe and the scalar field only depending on 
t, the FRW equations corresponding to (IHT!) and fl82!) are: 

= - ^H' + V{<P) + 24H'^^^^^ , (101) 

dt 



\ ) 



df^ 



■ deiUit)) odfl . ^ 

- 16HH " - im^ ,^ " . 102 

dt dt ^ ^ 

Combining the above equations, one gets 

dt^ dt dt 

.l^^_84f£!*(M)) , ,103) 

k2 dt\ a dt ) ' ^ ' 

which can be solved with respect to ^i(0(t)) as 

Combining fllOll) and fll04p . the following expression of V{(f){t)) may be found: 

VUit)) = ^Hitf - 3a(t)H(t)W(t) . (105) 

Due to a freedom of the redefinition of the scalar field we may identify t with 0. Then 
one considers the model where V{(j)) and ^i{(f)) can be expressed in terms of a single 
function g as 

\/(0) = -g' {<j>f - 3g' (0) e^(^)f/(0) 6(0) = tt / ^'/'i^7TT?f^('/'i) > 



n'"'^' .v.. v.-. ..v.-. "^(?'(0l)2 

f/(0) = ^ rc?0ie-^('^^V(0i) • (106) 

By choosing 1^(0) and C,i{4>) as fll06p . we can easily find the following solution for 
Eqs.dMD and (Hm : 

a = aoe^W {H = g\t)) . (107) 

Then we can reconstruct F(G')-gravity in the way very similar [351 EZ] to the scalar- 
Gauss-Bonnet theory in the previous sub-section. 

Although the above formulation is very similar to that in the scalar- Gauss-Bonnet 
theory, there could be some difference. In the scalar- Gauss- Bonnet gravity, since the 
scalar field has a kinetic term, the scalar field could propagate and there could be a 
possibility to generate extra force besides the Newton force. On the other hand, F{G)- 
gravity has no kinetic term for scalar field and extra force could not be generated. In 
fact, one can consider the perturbation around the deSitter background, writing the 
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metric as (^^j, = g(^o)f^u + hfj_^. Here the Riemann tensor in the deSitter background is 
given by 

R{0)fiiypa = Hq [g{0)^pg{0)ucj - g(0)tiag{0)i^p) ■ (108) 

The flat background corresponds to the hmit of i^o ^ 0- Fot simplicity, if we choose the 
gauge conditions g'^Q^h^^ = V^^-^h^^ = 0, we find from the equation of motion without 
energy-momentum tensor, 

= ^ V - 2iJoV) • (109) 

Since the contribution form the Gauss-Bonnet term does not appear except the length 
parameter I/Hq of the deSitter space, the only propagating mode should be graviton in 
the F(G)-gravity. This also shows that Newton law is effectively the same as in GR. 

Finally, let us note that it is not difficult to extend the above reconstruction scheme 
for modified Gauss-Bonnet gravity so that the (ideal fluid) matter may be naturally 
included into the formulation [37] . 

6. Discussion 

In summary, the reconstruction program is developed for the number of modified 
gravities: scalar-tensor theory, f{R), F{G) and scalar- Gauss-Bonnet gravity. It is 
explicitly demonstrated which versions of above theories may be reconstructed from 
the known universe expansion history. Specifically, it is shown that the cosmological 
sequence of matter dominance, decceleration-acceleration transition and acceleration era 
may always emerge as the cosmological solutions of such modified gravities. Moreover, it 
is explained that the (exact or approximated) ACDM cosmology may be also the solution 
of such gravity theories. Several examples of corresponding reconstruction for it as well 
as for oscillating universe where ACDM dark energy describes the one of the oscillation 
periods are given. In principle, it is not difficult to include into this reconstruction 
scheme also radiation dominated era and, perhaps, the inflationary epoch. This will be 
studied elsewhere. 

It is interesting that even if specific modified gravity which is well suited with 
acceleration epoch does not describe the matter dominance period, such a period may 
be included there as the solution at the price of the introduction of compensating dark 
energy. The corresponding example is worked out for f{R) model of ref. [23]. It may be 
also extended for scalar- Gauss-Bonnet or F{G) gravity. 

It is quite remarkable that modified gravity may not only naturally describe dark 
energy (unlike to conventional GR) but also it can be reconstructed using the realistic 
universe expansion history. It is also promising that it may be fitted with observational 
data. As well it may be compatible with Solar System tests. Of course, some 
controversial results exist here what is not strange due to fact that modified gravity 
is seriously considered as the alternative cosmological theory only last several years 
while GR played this role almost one century. One should also bear in mind that most 
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of cosmological data were obtained (at least, up to some extent) using GR foundation. 
Hence, very serious reconsideration (as well as more precise and complete observational 
data) are requested in order to have the answer to the fundamental question: what is 
the gravitation theory which governs the expansion of our universe? 
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